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Abstract 
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The profile of a cylindrical shell in the post-buckling regime of deformation is calculated exactly, for any load and axial deformation, 
and is shown to be a Jacobi elliptic sine function. Cylindrical symmetry is assumed from the beginning. Results shed light on why 
f^**) old tires are so effective as dock and tugboat bumpers from long ago to nowadays, proving to be safe for berthing even the largest 
freighters. As no technical analysis of this capacity has been published yet, it seems opportune to apply the mathematical solution 
Q to explain this long lasting enigma. The reaction force and stored energy of an axially compressed cylindrical shell, as a tire tread 
<D is, exhibit ideal behaviour for a safety shock energy absorber. Most energy absorption takes place after buckling, which bends the 
cylindrical mantle and increases its circumference, inducing strong tensile forces that contribute to stabilize the buckled structure. 
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1. Introduction 

There is a great deal of work done on the critical conditions 
for the buckling of narrow or thin walled solids of a variety of 
shapes, and on the subsequent evolution of the collapsing ele- 
ment [ 1 1 . Buckling is generally the controlling failure mode of 
thin shell structures and long narrow structural elements, hence 
the study of the elastic instabilities of shells is a most important 
subject of materials mechanics and engineering design. There 
is an extensive list of books, review articles and conference pro- 
ceedings providing a wealth of information on the strength, me- 
chanical stability and buckling behavior of thin shell structures 
0. However, either explicitly or not, most of these studies 
have associated the idea of buckling as a catastrophic collapse 
of the critically strained thin walled or narrow element. Also, 
the mathematical treatments are often intended to be as gen- 
eral as possible |3j [4] [5] 17), which precludes the obtention 
of closed-form practical formulas of use in specific practical 
situations, and able of being tabulated in manuals (TJ. 

This communication deals with a mechanical model for a 
very particular technical application in which buckling is not an 
undesired effect, but a highly beneficial one. The use of old tires 
as dock bumpers has been an old tradition in most harbours. Al- 
though they are now being replaced by specially fabricated ma- 
rine fenders, ship captains appreciate the behaviour of old tires 
in docking because, for instance, they provide a characteristic 
smooth absorption of the ship kinetic energy with total absence 
of spring back at the end of the operation. The question of why 
old tires are still being preferred as dock protections, even for 
the docking of large cargo vessels, is an enigma whose precise 
solution is the subject of what follows. 
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The reaction force and energy absorption effect of an axially 
compressed tire is attributed entirely to the elastic deformation 
of the tread, which is modeled here as a cylindrical shell whose 
radius R is greater than its axial length Lq, both magnitudes 
measured in the unstrained condition, as shown schematically 
in Fig.[T] In opposition to small tires, whose treads are often ini- 
tially slightly barreled, big ones are essentially cylindrical when 
unloaded. Though small, the shell thickness e is assumed large 
enough to preclude the formation of patterns of complex geom- 
etry upon deformation. The applied opposed forces of strength 
P are parallel to the symmetry z-axis and are uniformly dis- 
tributed along the edges of the cylindrical mantle. It is also 
assumed that the edges of the elastic cylinder are not deformed 
because they rest on rigid plane surfaces and the friction forces 
are strong enough to prevent sliding. Hooke's law is assumed 
to hold in a first approach. Although rubber elastic properties 
exhibit deviations from linear response, the main conclusions 
will not be much affected by this. 

Instead of dealing with the model in the most general way, 
facing all the complex deformation modes a cylindrical shell 
may display |3] 01 [5] |6] 0, the mathematical approach is kept 
as simple as possible by allowing only barrel shaped deforma- 
tions. The equilibrium equations of the model incorporate from 
the very beginning just the deformation modes one can observe 
in an axially compressed tire, and discards any other because 
the main objective is to understand the ability of tires as shock 
energy absorbers. 

The mechanical analysis shows that, in a short first stage, the 
cylindrical shell undergoes a uniform compressive strain along 
the z-axis, conserving strictly its cylindrical shape and oppos- 
ing a reaction force proportional to the strain up to a maximal 
load Pb- The maximum strain reached in this first deformation 
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regime is indicated in Fig. [TJa) by a discontinuous line. Com- 
pression beyond this limit makes the shell to start buckling, and 
the cylindrical mantle progressively acquires a barrel shape, as 
in Fig.[TJb). In this second regime the reaction force diminishes 
mono tonic ally with strain, but conserving always a significant 
strength. The development of the barrel profile involves stretch- 
ing of the tread circumference, which demands strong tensile 
forces that contribute to stabilize the buckled structure. As a 
shock energy absorber, the system exhibits the advantages of 
having a maximum reaction force, which can be designed to 
warrant no damage to the colliding bodies, together with a large 
energy absorption capability, particularly in the buckled strain 
regime. 



z 




Figure 1 : Model for the tread of a pneumatic tire as a cylindrical shell, (a) 
Unstrained condition (solid lines) and critical axial strain for buckling (discon- 
tinuous lines), (b) Buckled strain regime by the action of compressive axial 
forces homogeneously distributed on the shell edges. 
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Figure 2: The strained cylindrical shell, showing schematically an elementary 
sector defined by the coordinates r, <f> and z, (0 is implicit), and by the variations 
Aip and Az. 



schematically represented in Fig. [2] The method is preferred 
over simply writing the equations of the theory of elasticity and 
applying boundary conditions because is simpler and provides 
better physical insight. The advantages of the procedure are 
discussed in the final subsection of this section. 

Fig. [3] displays a diagram of the material element and the 
forces applied on it. Forces in the plane <p = constant, contain- 
ing the z-axis, are compressive forces of strength F c (z + Az) 
and F c (z). The forces in the plane z - constant, normal to the 
z-axis, are tensile forces because stretch the shell contour in 
its plane, and their strength is denoted F,p(z). Solutions with 
cylindrical symmetry were implicitly assumed because no de- 
pendence of the forces on <p was considered. 




' F c {z) 

Figure 3: The elementary sector of Fig. [2] and the compressive and tensile 
forces F c and exerted on it. 



2. Buckled strain regime of the cylindrical shell of finite 
thickness 

2.1. Equilibrium of the internal forces 

Fig. [2] shows the cylindrical shell in the strain regime in 
which it adopts a barrel shaped profile by effect of the axial 
load P (equivalent to the reaction force at equilibrium, or quasi- 
equilibrium). The cylindrical coordinate system has its z-axis 
along the main symmetry axis and the origin O at the center of 
the deformed shell, which extends from z = -1/2 to z = 1/2. 
Distance r is the radius of the shell for a given z, so that func- 
tion r(z) - R determines the profile of the strained shell. The 
azimuthal coordinate <p is the usual one in cylindrical coordi- 
nates, and is not shown in Fig.[2]for the sake of clarity. 

The equilibrium equations are derived from the analysis 
of the forces operating on an elementary sector of the shell, 



Fig. Qa) is a projection on the plane z = constant of the 
system of Fig. [3] showing the forces operating in this plane, and 
Fig. |4|b) shows the vector composition of them. Hence, the 
tensile forces give a sum of strength 

2F $ {z) sin(A0/2) ^ F (z)A0 if A<f> * 0, (1) 

contained in the plane z = constant and pointing towards the 
central z-axis. 

Forces in the plane <p = constant are shown in Fig. [5] They 
are the two compresive forces of strength F c (z + Az) and F c (z) 
operating in the two opposite edges, and the third one is the 
resultant of the tensile forces in the plane z = constant, whose 
strength is F ( f ) (z)A(f). The former two forces are tangent to the 
curve r(z)-R determining the shell profile in the corresponding 
application points, and subtend angles 9 + A8 and 8 with the 
z-axis, respectively. 
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Figure 4: (a) The graphic scheme of Fig.[3]projected on the plane z. = constant, 
(b) Force diagram in the plane z = constant. 
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Figure 5: Forces in the plane (p = constant. Compressive forces F c (z + Az) and 
F c (z) are tangent to the curve r(z) - R representing the shell profile, and subtend 
angles 8 + A8 and 8 with the z-axis, respectively. The force F*(z)A(fi is obtained 
in Fig.[4]from the tensile forces in the plane z. = constant. The discontinuous 
line represents the shell profile. 



Equilibrium demands 

- F c (z + Az) cos(0 + A0) + F c (z) cos = 0, (2) 

F c (z + Az) sin(0 + A0) - F c (z) sin = A0. (3) 

From Eq. ([2]) one can infer that the z-component of the forces 
does not depend on z, and hence they can be identified with the 
external force PA(f>/(2n) applied on the edge of the shell fringe 
defined by A<p. The argument is graphically shown in Fig.[6]and 
yields 
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Figure 6: Equilibrium of a shell fringe of finite size in the z-direction. The 
z-component F c {z) cos 8 of the force is PA<f>/(2n) for any z in [-1/2, 1/2]. 



tan(0 + A0)-tan0 = 
Force F$ can be expressed as 

F$ = cr^eAs, 



2nF& 



(5) 



(6) 



where cr^ is the normal stress in the azimuthal direction, e is 
the shell thickness and As the length of the shell element shown 
in Fig. [5] and whose value is 



As = x / 1 i-l^l A". 



(7) 



Replacing tan0 = dr/dz and combining Eqs. (|6]l and (|7]i 
one obtains that, in the limit Az — > 0, 



(8) 

and cr 77 sat- 




On the other hand, the normal stresses cr rr , 
isfy Hooke's law 



r-R 1 

£00 = — — = -(O"00 - VO - V(T ZZ ), 

K h 



(9) 



where E and v are the Young coefficient and Poisson ratio of 
the material the shell is made of. Substituting 



2nre 



cos 0, 



(10) 



cos0= [1 + (dr/dz) 2 T l/2 and 

F$A<p 

o~ rr — = 

rA^Az 

in Eq. d9]l, and solving for <r^, it is obtained that 
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F c (z + Az) cos(0 + A0) = F c (z) cos = — Ad>. (4) 

2n 



2.2. Differential equation for the shell profile function r(z) — R 
Recalling Eqs. (HI) one can divide the first term in the right 
hand side of Eq. Oil by F c (z + Az) cos(0 + A0), the second one 
by F c (z) cos 0, and the right hand side by PA4>I(2tt). It gives 
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(12) 



Combining Eqs. ( 12 1 and ((8j) one finally obtains 



dz 2 



2neE 
~RP 



(r-R)Jl + 
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R 



1 + v- 



(13) 



The differential equation (13i gives the profile r(z) of the 



cylindrical shell deformed by the axial load P when the proper 
physical conditions are imposed to the solutions. The most ev- 
ident conclusion one can infer from Eq. ( p"3j ) is that the load P 
cannot be null. Hence it corresponds to a deformation regime 
occurring for high enough loads. 

The unloaded cylindrical shell has a length L in the axial di- 
rection, which reduces to a smaller length / by the applied load 
P. The length reduction originates from the elastic compressive 
strain of the deformed mantle and by the shape deformation 
itself projected in the z-direction. If the elastic strain can be ne- 
glected when compared with the effect of the shell deformation, 
the length of the curved path adopted by the deformed shell will 
conserve its original value L. Hence the physical solutions of 
Eq. ( 13 1 must satisfy 



drY 
d~zj 



(14) 



2.3. Present theoretical approach versus theory of elasticity 

The equation for the shell profile was derived in the preced- 
ing subsections from analysing the equilibrium of the forces 
exerted on a representative elementary sector of the material 
medium, introducing from the beginning the symmetry con- 
straints. For our present purposes, this procedure is much sim- 
pler and practical than the more standard approach of readily 
introducing the system boundary conditions into the general 
formalism of the theory of elasticity. This is because the lat- 
ter approach oblige us to face complex problems we are not 
interested in. In the theory of elasticity the free surfaces of 
the shell are defined as surfaces where the stresses vanish, and 
the theory furnishes the equations for calculating the detailed 
structure of the stresses occurring in between. A recent paper 
by Zozulya and Zhang [8 1 provides a good example of this kind 
of precise calculation of the deformation of cylindrical shells, 
together with a detailed account of the stress fields inside the 
material. The cost paid for such a complete solution is the in- 
troduction of numerical methods from the very start. 

However, it is not necessary here to know about the precise 
profiles of the outer and inner surfaces of the shell, or how vary 
the stresses inside. The forces F c {z) and F^,(z) in Fig. [3] account 
for the integrated effect of these stresses, and their equilibrium 
condition proves to be enough for determining the equation for 
the mean shell profile, which has the advantage of being closed- 
form. 



3. The profile equation in the buckling regime of deforma- 
tion and its solution 



(r = constant), 



(15) 



R 2nre 

which can be interpreted as the Poisson effect on the shell circu- 
lar perimeter 2nr accompanying the homogeneous axial strain 
produced by the applied compressive stress P/(2nre). This 
uniform solution, which is expected to be stable up to a criti- 
cal load Pg, preserves the cylindrical shape and produces only 
small geometrical variations, even for big loads, because E 
is usually very large (0.01-500 GPa). The buckled non triv- 
ial solutions of Eq. ( p"3} involve much larger deformations. 
In the post-buckling regime r assumes values in the interval 
R < r < R + L/2 when I varies from / = L to I = 0. The second 



term v/r in the numerator of the right hand side of Eq. ( 13 i is 
comparable with the first one only when \r - R\/R <k 1. As 
long as r departs from R beyond the range of the purely elastic 
distortions, the term v/r becomes negligibly small when com- 
pared with the one proportional to E. Hence it is advisable to 
distinguish between elastic strains and geometric changes and 
write 



d 2 r 
dz 2 
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(buckling regime). (16) 



Ec. ( [T6| > also assumes a thin enough shell to neglect the second 
term in the denominator of the right hand side of Ec. ( 13 1. Al- 
though this non-linear equation is not listed in the specialized 
treatises on elliptic functions and integrals ||9] [10), it will be 
shown next that its exact solution is a Jacobi elliptic sine func- 
tion, which holds for any load and deformation state, including 
the limit in which the shell has been completely flattened by the 
applied force. 
Defining 



2neE 

y(z) = - R, p — , R = constant, 



PR 



Eq. ( 16 1 reads 



and the substitution 



turns it into 



= -py V 1 + y' 2 > 



y(z) = — p u{ yp z) 



u" - —u Vl + u' 2 . 



(17) 



(18) 



(19) 



(20) 



Multiplying both sides of this equation by u' yields the inte- 
grate form 



u u 



Vl + u' 2 
which can be solved to give 



= — uu , 



The trivial solution r = constant reduces Eq. ( 13 i to 



Vl + U' 2 = --U 2 + -Kn + 1, 

2 2 



(21) 



(22) 



where uq = u(0). Because of the symmetry with respect to 
the origin, u( yfp z) must have its maximum at z = and the 
integration constant was chosen so that u'(Q) = 0. Denoting 
now 



Z=y[pz, v(0 



Uq 



(uq/2) 2 
1 + («o/2) 2 



Eq. ( 22 1 can be rewritten as 



dv 



or 



Vl -k 2 



dv 
~dl 



V(l - v 2 )(l -k 2 v 2 ). 



(23) 



(24) 



(25) 



Inverting this equation and integrating with respect to v one has 



1 



Vl -k 2 



,»=r 

Jo 



dv 



V(l -v 2 )(l -k 2 v 2 ) 



= F(v,k), (26) 



where F(v, k) is the incomplete elliptic integral of the first kind 
with modulus k, (0 < k < 1) |0[TO), and £\ is an integration 
constant. The same symbol was used for the integration vari- 
able and the upper integration limit to simplify the notation. 

The inverse function of the incomplete elliptic integral 
F(v, k) is known as the Jacobi elliptic sine function sn. Hence 



Eq. ( 26 1 is equivalent to 



sn 



(27) 



Function sn (x, k) takes values in the interval [-1,1] and is pe- 
riodic in x with period AK{k), where K(k) = F(l,k) is the 
complete elliptic integral of the first kind. Also, sn (0, k) = 
sn (2K, k) — and sn (K, k) — 1 . The Jacobi sine function is 
symmetric with respect to x = K and has a maximum there. 

Therefore, the solution satisfying the boundary conditions 
y(±l/2) = is such that 



1 -k 2 



2K(k) 



and 



Zl = 



(28) 



Taking -Jp and z\ from these equations, and recalling the third 
of Eqs. (23 1, which gives mq/2 = k/ Vl - k 2 , one has that the 
profile of the buckled cylindrical shell is given by 



y(z) = r(z) - R 



(1 -k 2 )K(k) 



sn 



2K{k)( /< 
— l Z+ 2, 



(29) 



Combining the first Eq. (|28) with the definition ( 17 1 of the 



be shown that k 2 is essentially the relative axial deformation 
e — (L— 1)1 L of the cylindrical shell. The elliptic integral K(k) 
can be calculated quite easily from the defining integral or the 
series 



(31) 



o sin(Tix) 
— sn(2/fx) 




Figure 7: The Jacobi elliptic sine function (solid lines) for three values of 
the modulus k, compared with the trigonometric sine function (circles). For 
k < 0.7 the two functions are very similar. At k > 0.7 the elliptic function 
is sensibly broader and goes to a square wave in the limit k = 1 . Physically, 
k = 1 corresponds to the situation in which the cylindrical mantle has been 
completely flattened by the applied force P and / = 0. 



Fig. [7] shows the Jacobi sine function sn (2Kx, k) for three 
values of k. The graphs makes apparent that the Jacobi elliptic 
function goes from sin(7rx) to a square wave when the modulus 
k goes from to 1 . The latter situation corresponds to the final 
collapse of the cylindrical mantle, when I = 0. 



3.1. The equation for the modulus k 
When R = constant one has r'(z) ■ 

2x1/2 



y'(z) = u'( t/pz), hence 



constant p, arelation between the applied force P and k follows ( l + r ' Z )' /Z can be replaced by (1 + M ' 2 )^ 2 in the general con- 



neEl 2 



, , -■ (30) 
2R{\ -k 2 )K 2 (k) V ; 

Therefore, it rests just to determine the meaning of k to have a 
complete solution of our problem. In the next subsection it will 



dition expressed by Eq. 
Eq. ((22) 

~//2 



14 1. Making this and then substituting 



rl/2 i W/2 _ : 

L- I Vl + u' 2 dz = -- I u 2 (yfpz)dz+ -u 2 Q l + l. 

J -1/2 2 J_i/2 2 



(32) 
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Recalling now 



2K( l\ 
-( Z+ 2)> k 



2k 



Uq 



Vl + k 2 



(33) 



Eq. ( 32 1 can be written as 



L-l 



2k 2 
l-k 2 



sn 2 (£,k)d£ 



(34) 



The integral appearing in this equation has been solved in terms 
of the elliptic functions and one has the mathematical identity 

m 



s 



! ^ = i[«-E(ama)], (35) 

where E(<p,k), < (f> < n/2, is the second kind incom- 
plete elliptic integral with modulus k. In the usual notation 
of the theory of elliptic integrals the amplitude am means 
am £ = arcsin(sn 0. Care must be taken in replacing prop- 



erly the integration limits in the indefinite integral ( 35 i because 
< am £ < n/2, and am K = n/2. Therefore, 2K is outside the 
domain in which the amplitude function is defined. To over- 
come this difficulty one can take advantage of the symmetry of 
sn (x, k) with respect to x — K and write 



which determines exactly the modulus A; as a function of only 
the relative axial compression s = (L — 1)1 L. The function k = 
/(e) does not depend on the characteristics of the cylindrical 
shell and is shown in Fig. [8] The amazingly simple approximate 
rule 



k 2 * e 



(42) 



holds over the whole range of e with a 5% maximum error. 
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Figure 8: The squared modulus k 2 and the relative axial compression e 
(L - t)/L, as given by Eq. |5T). 



sn 2 <T^ = 2 sn 2 ^ =-r[#- E(n/2,k)] (36) 
Jo Jo k 2 

replacing am K = nj2. E(n/2, k) = E(k) is the complete elliptic 
integral of the second kind, modulus k, defined by the integral 
13 



E(k) 



f 



or the series 



2 sin (p dip, 



2 [ 2"nl I 2n + 1 



(37) 



(38) 



The term D(k) = (l/k 2 )[K(k) - E(k)] can be evaluated either by 



combining the power series ( 3 1 1 and ( 38 1 as 



D(k) = 



K(k) - E(k) = n y 2{n+\) l (2» + l)H y 2n 
2 2n+ 1 \2" +1 («+ 1)!/ 



11 \2 



k 2 



or solving the integral 



D(k) 



■"12 cin 2 



sin (p dip 



° ^|^k~ 2 



(39) 



(40) 



sin 



which follows directly from the definitions of K(k) and E(k). 
Eq. ( |34| then becomes 

£ fl-^l s= L -^l, (41) 



1 - s l-k 2 \ K(k) 



L ' 



4. Critical load and reaction force 

The deformation regime studied in the preceding sections 
takes place for a constant non-vanishing axial load P, given 
by Eq. (30 1. This is because it was assumed from the very 
beginning that loading curves the cylindrical shell profile and 
Eqs. (29 1, ( 30 » and (41 1 hold at buckling and in the post- 
buckling regime of deformation. The variable e measures only 
the deformation induced relative displacement of the edges 
along the axial direction, and the much smaller elastic distor- 
tions have been neglected. The force 



neEL 2 (\ - ef 
2R(\ -k 2 )K 2 (k) 



(43) 



has an absolute maximum Pb at s = k — 0. Recalling K(0) = 
n/2, the maximal reaction force reads 



and one can write 



Pie) 
Pb 



2eEL 2 
nR 



(1-s) 2 



(44) 



4 (1 -k 2 )K 2 (k) 
Fig. [9] shows P/Pb, which according to Eqs. ( |45] l and (41 



(45) 



an universal function of s. Fig. 10 displays the energy stored 



by the shell in the post-buckling deformation regime, obtained 



integrating Eq. (45 1 with respect to s 
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Figure 9: Reaction force P relative to the maximum Pg as a function of axial 
compression e = (L — l)/L , as given by Eq. \45). 




0.2 0.4 0.6 0.8 
Relative axial compression £ 



1.0 



Figure 10: Energy absorbed by the cylindrical shell in the post-buckling de- 
formation regime, in units of PgL and for any relative axial compression e. 



5. Homogeneous and buckled deformation regimes 

The problem we are interested in is a deformation in which 
the cylindrical shell is gradually compressed along its main 
symmetry axis. In a first stage the compression is elastic and 
conserves the cylindrical shape. The relative displacement x of 
the edges is given by Hookes's law 



7lR 



nR 



nR 



\ y/2Lo) V2L \\ V2L 



+ 2, 



(47) 



which is the equation for the critical elastic strain x b /Lq. The 
critical load is 



P B = InReE 



1 + 



nR 



nR 



nR 



( V2Z.J V2I,, VI \'2L, t 



(48) 



6. Conclusions relative to the marine application 

To illustrate why vehicle old tires are functional as dock 
dampers in the light of the preceding equations, consider a 
27"x 47" tire, originally designed for heavy mining trucks. The 
dimensions of the tread are Lq = 0.75 m, R = 1 .25 m and 
e = 0.10 m. Young's effective modulus for tire rubber with 
added reinforcements can be estimated as £ = 10 MPa. The 
purpose of this calculation is fundamentally illustrative, to ac- 
quire a feeling of the magnitudes involved. Replacing these 



constants in Eqs. (47 1 and (48 i one has that 



x B = 2.55cm P B = 267 kN = 27.3 ton. 



(49) 
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Figure 1 1 : Axial reaction force P exerted by a tire with tread dimensions 
Lq = 0.75 m, R = 1.25 m and e = 0.10 m, and effective Young's modulus 
E = 10 MPa, as a function of the axial compression x. The broken line shows 
the behaviour of a spring or any other strictly elastic medium chosen to yield a 
maximum force comparable with the tire. 



P X 

= E-, 

2nRe L 



(46) 



where Lq is the original length of the unstrained cylinder. Once 
the elastic distortion x takes a critical value xb, such that the 
applied force P reaches the threshold strength P B , the cylinder 
buckles and Eq. (43i starts holding, instead of Eq. (46 1. At 
the critical deformation x B the surface just starts to acquire the 
barrel shape, the two deformation regimes coexist and Eqs. (43 1 
and ( |46*1 ) are both valid for L - Lq - x B and e = (notice that L 
stands for the the shell length when buckling is just initiated). 
Thus the axial deformation x B at which the buckling regime 
sets in can be obtained from combining Eqs 
eliminate P B . This yields 



44 1 and (46 1 to 



Thus, when the tire is progresively compressed, the reac- 
tion force grows linearly with strain up to the critical strength 
P B — 267 kN at x B — 2.55 cm, with the tread conserving its 
cylindrical shape. As compression continues, the tread curves 
and the radius turns dependent on z according to Eq. ( |29) . 
Fig. 1 1 shows the behaviour of the reaction P of the strained tire 
in the two regimes, the linear and the buckled one (solid line). 
In the latter regime force P obeys Eq. (45 1 with L — Lq - x B . 
The broken line is the force exerted by a strictly elastic system, 
like a spring matrix, introduced here for the sake of comparison. 

The stopping power of the tire when used as a dock bumper 
device is given by the area W{x) under the curve of the reaction 
force. The solid line in Fig. 12 shows the curve of the energy 
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Figure 12: The solid curve represents the energy absorbed by the tire of Fig. fTT] 
when compressed in a distance x. The broken curve is the same magnitude for 
the strictly elastic system. The smaller capacity of the latter for compressions 
in the range between and 40 cm is apparent. 



stored W(x) by the tire as a function of the compressive dis- 
placement x. The broken line displays the same magnitude for 
the purely elastic shock energy absorber. 



Figs. 11 and 12 show the features that make tires good dock 
bumpers: 

(a) The reaction force has a well defined maximum which 
ensures that the energy absorption device will not damage the 
colliding structures. 

(b) Once the maximum reaction is reached, the reaction P de- 
creases smoothly with the displacement x, and energy storage 
W(x) keeps high. This is the physical magnitude determining 
the stopping power of the device. 

(c) Function W(x) increases monotonically and rapidly with x, 
particularly for small values of x. This means that the stop- 
ping action of the axially compressed tire rapidly increases and 
keeps high at any strain, even for small ones. The superiority 
of the cylindrical shell (solid line) over a strictly elastic system 



(broken line) is shown quite dramatically in Fig. 12 The latter 



reaches an acceptable stopping ability at the cost of increasing 
the reaction force to a dangerous level. 

(d) As x takes relatively large values, the reaction force P of the 
cylindrical bumper decreases to a fraction of its initial strength, 
reducing this way its ability to produce rebound. The elastic 
bumper produces the opposite effect. 

Think for instance of a ship of mass m = 55, 500 ton being 
pushed by the tugs to approach normally the dock at a speed 
v = 0.24 knot = 0.1234m/s (a rather high speed). The kinetic 
energy of the ship is then 423 kJ. Assuming a safe tire com- 



pression of 0.4 m, the graph in Fig. 12 shows that each tire can 
store W(0.40m) = 72.5 kJ, and then a minimum of six units 
would be necessary for safely docking the vessel. On the other 
hand, the ship side is specified to support a maximal load of 
147kN/m 2 . The peak pressure the tires can exert on the ship 
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side is P B /(nR 2 ) 
fied maximum. 



54.4 kN/m , which is close to 1/3 the speci- 
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